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Abstract 

A collapsing spherical D2-brane carrying magnetic flux can be described 
in the region of small radius in a dual zero-brane picture using Tseytlin's 
proposal for a non-Abelian Dirac-Born-Infeld action for N DO-branes. A 
standard large N approximation of the DO-brane action, familiar from the 
Myers dielectric effect, gives a time evolution which agrees with the Abelian 
D2-brane Born-Infeld equations which describe a D2-brane collapsing to zero 
size. The first 1/N correction from the symmetrised trace prescription in the 
zero-brane action leads to a class of classical solutions where the minimum 
radius of a collapsing D2-brane is lifted away from zero. We discuss the 
validity of this approximation to the zero-brane action in the region of the 
minimum, and explore higher order 1/N corrections as well as an exact finite 
N example. The 1/N corrected Lagrangians and the finite N example have 
an effective mass squared which becomes negative in some regions of phase 
space. We discuss the physics of this tachyonic behaviour. 
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1 Introduction 



There has been a lot of work recently on time dependence in string theory 
[H 12 El HI 13 El H|- In this paper we will study a simple time dependent 
system of a spherical D2-brane carrying magnetic flux. 

D2-branes in type IIA string theory with magnetic flux on the world- 
volume carry zero-brane charge [HI E] ■ We consider such spherical D2-branes 
and their time evolution as described by the 2+1 dimensional Born-Infeld 
action. A closely related problem is the time dependence of a spherical M2- 
brane of M-theory. When considered in the context of Matrix Theory ^21 EI] > 
we are led to add momentum along the eleventh direction, which amounts 
in the IIA picture to having D2-branes with magnetic flux. Other closely 
related systems include the microscopic description of giant gravitons [12] 
and the D3-brane bion where the world-volume of the D3 stretches out into 
a Dl-funnel described by an excitation of the scalars in the D3-world-volume 
|13[ HI] . Magnetic fluxes in the D3-worldvolume allow the existence of a 
BPS solution describing such funnels. The dual description in terms of the 
Dl-worldvolume has been studied in detail in [IS], one general lesson being 
that the Dl and D3 descriptions agree at large N. Whereas the D1-D3 funnel 
system involves a spherically symmetric solution of the Dl-world- volume with 
the radius R(cr) of the spherical cross-section of the funnel having a non- 
trivial dependence on the spatial coordinate a, the D0-D2 system has a time 
dependent radius R(t) of the spherical D2-brane. 

In a similar spirit to HJ>] we may look for the collapsing D2-brane in terms 
of the DO-branes. The dualities between descriptions of the same physics from 
two points of view, of a lower dimensional brane and a higher dimensional 
brane, can thus be explored in a time-dependent context. The DO-brane 
effective world-volume Yang Mills theory neglects stringy excitations. This is 
a valid approximation when the separation of the zero branes is less than the 
string length (branes as short distance probes in string theory were studied 
in detail in J7]). For a spherical system of N DO-branes this means that 
the radius R of the sphere obeys R < l s VN. The semiclassical world- volume 
D2-brane action is expected to be valid for R > l s . There is a large overlap of 
regimes of validity at large N. Indeed we will find in Section 2 that the non- 
Abelian Born-Infeld type action of the DO-branes (written down in and 
developed to describe the dielectric effect in [!§]) gives at large N the same 
equation for the time evolution of R as the one obtained from the D2-brane 
action. 

In Section 3, we consider 1/N corrections to the equation for R. This 
requires some combinatorics of symmetrised traces. We give two methods 
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for computing these corrections. The first is based on evaluation of a sym- 
metrised string of su(2) generators on the highest weight state of an irre- 
ducible representation. The second method maps the problem of expressing 
the symmetrised generators in terms of powers of the Casimir, into a combi- 
natorics of chord diagrams of the kind that appear in knot theory. 

In Sections 4 and 5 we study the leading 1/N correction. We consider 
the energy as a function of (r,s), the dimensionless position and velocity 
variables. For convenience many formulae are expressed in terms of the 
variables (7,?/) where 7 = (1 — s 2 )~2 is the standard Lorentz factor of 
special relativity, and U = (1 + r 4 )a. At leading order all the formulae 
for Lagrangian, energy, momentum, etc, are the standard ones of special 
relativity, with U playing the role of a position dependent mass. The 1/N 
corrections give an interesting modification of these standard formulae. The 
physical consequence of the correction is that while a spherical brane starting 
at rest at a large radius collapses to zero radius for a range of radii, it cannot 
do so if the initial radius is too large. These large collapsing membranes 
bounce from a minimal radius, which depends on the initial radius. This is a 
somewhat exotic bounce where the velocity r is discontinuous, as in reflection 
at a hard wall, and in addition the fate of the membrane after the bounce is 
not uniquely determined, but rather there are two possible outcomes. These 
conclusions follow from an analysis of the energy function in (r, s) space and 
the existence of an extremum where §f = and s ^ 0. We briefly discuss 
quantum mechanics near this extremum. 

In Section 6 we consider higher order corrections in the 1/N expansion, 
explaining which properties of the solutions are preserved as we include these. 
In Section 7 we consider the exact finite N evaluation of the symmetrised 
trace for the case of spin 1/2. We find explicit forms of the energy function in 
terms of hypergeometric functions and discuss the dependence upon position 
and velocity. In Section 8 we discuss regimes of validity and the effects of 
higher order terms in general, paying attention to the behaviour of the accel- 
eration and effective mass, and the appearance of tachyonic modes. Finally, 
in Section 9 we present a summary and outlook. 

2 Time dependent solutions: Born-Infeld 

Consider the action describing N DO-branes 

So = -T / dt Str^j (I N - \*d t &QT j 1 d t &) Vdet Q« (1) 
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where , i = 1 . . 9 are a set of N x N matrices and Q ij = I N S ij + T 
is the tension, A = lua' = 2nl 2 and Str denotes the completely symmetrised 
trace over iV-dimensional indices. Amongst the solutions to the equations of 
motion obtained from ([I]) are those describing t-dependent fuzzy 2-spheres. 
These are obtained from the ansatz = R(t)oti, i — 1, 2, 3; <£> m = 0, m — 4. .9, 
where a% are the generators of su(2) in an iV-dimensional representation, with 

[ai,aj] = 2it ijk a k (2) 

The physical radius of the fuzzy S 2 which we denote by R(t) is related to R 
via 2 

R 2 = ^tr(® 2 ) = \ 2 CR 2 (3) 

with C = (N 2 — 1) is the value of the second order Casimir invariant. 

Substituting this ansatz into the action (JTJ) and using the approximate 
relation (valid in the leading large iV limit) Str(diai) n = NC n , we get 

S = -T N J dtyj (1 - {\ 2 C){d t R) 2 ) (1 + A\ 2 CR A ) (4) 



The conserved energy is E = « / 1+4X2 c ^ 4 an d the equations of motion are 



Integrating this equation with the initial condition R = R and R = at 
t — 0, we obtain 

where we work in physical units. This equation, studied in the context of 
time-dependent membranes in [20], has solution 



R(t) = R Cn(tV2/R , 1/V2) (7) 

where Cn(u,k) is a Jacobi elliptic function and R% = R ° + ^ 2 C ^ 4 ■ The initial 

conditions are satisfied due to the properties Cn(0, fc) = 1 and dCn j^' k ^ = 
—Sn(u, k)dn(u, k) with Sn(0, k) = 0. 

The functions Cn(u, k) have the property that they are monotonically 
decreasing with zeros at the special value u = K(k) where K is a complete 
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elliptic integral of the first kind. Thus our D2-brane solution describes a 
spherical collapse to a point after a time 

^§1 = K(l/V2) (8) 

ie the collapse time is = - (yi% + A 2 C/4) K(j^). Hence we see that 

increases as we increase the size of R (assuming that R$ » A 2 C/4). 
Solutions of a similar form have been discussed recently, in the context of 
3-form flux backgrounds, in [To"] 

Now consider the DBI action for D2-branes moving in flat space but 
including Abelian world- volume gauge fields - 

S 2 = -T 2 J d 3 ^-det(G MN d a X M d b X N + XF ab ) (9) 

where in © £ a , a = 0,1,2 are worldvolume coordinates of the D2-brane, 
Gmn is the flat space metric in D + 1 dimensions, M = . . . D, F ab is 
the Abelian gauge field strength and T 2 is the brane tension. We choose a 
static gauge where £° = t = X°, £ x = X 1 , £ 2 = X 2 . We want to consider time 
dependent solutions of a spherical D2-brane and so we shall choose world vol- 
ume coordinates £ = #,£ 2 = 0, and embedding X 1 = R(t) sin#sin0,X 2 = 
R(t) sin 9 cos <ft, X 3 = R(t) cos <j), where R(t) is the radius of the spherical D2, 
in physical units. 

In order to correctly reproduce the dynamics of the time dependent fuzzy 
spheres considered elsewhere, we shall take the field strength F ab to define n 
units of magnetic flux through the world volume, ie 

dL ab F ab = 2vm (10) 

where dT, ab is the infinitesimal world- volume surface element. Since only the 
F12 components are non-zero, F 12 = ns ™ ■ Its easy to see that the action for 
the D2-brane reduces to 



S 2 = -2nn\T 2 J dtj (l - i? 2 ) (l + 4-^) (11) 
The equations of motion for this sytem can be written in the form 

dt\ l-R* 1 
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Comparing the above equation to that describing a system of N DO-branes, 
discussed earlier, we see that in the leading N approximation the equations 
coincide if we identify N = n. 

It is useful to introduce dimensionless variables (r, s), where s — 4^, 
such that the Lagrangian is + r 4 )(l - s 2 ). The variable T is a re-scaled 
time t. The relations between the dimensionless variables and the original 
variables in (jlj) are 

r 4 = (AX 2 C)R 4 
s 2 = (X 2 C)R 

T = V2(\ 2 C)- l/ H (13) 

Using the relation between R and the physical R in (J3J) we also see that 

dr dR , , 

so that the relativistic barrier is just s = 1. 



3 Corrections from the symmetrised trace 

Consider the ansatz $j = Rcti, introduced in the DO-brane theory in the 
previous section. Substituting this into the action (JTJ) yields 

* = - r ./W(l + 4XW>)(l-AV«>) (15) 

with both % and j summed from 1 to 3. The derivation of ()15|) uses properties 
of the symmetrised trace. For the analogous discussion in the Dl-brane case, 
see The equations of motion following from (|THj) are 



Str- 
btr dt' 



l + Wa^R* Q (lg) 



\ l-\ 2 aiaiR 2 

In the large limit, we may replace a^aj by C(N), and we obtain (j^J). 
However, we wish to use more precise expressions here. To consider this, 
suppose that one wishes to evaluate the symmetrised trace of some function 
/(ajttj) of the summed product ajCtj. If we expand 

oo 

f(x) = Y,f^ n (17) 

n=0 
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for some constants f n , then 

oo 

StrffaiOti) = 22 fnStr(aiOii) n (18) 

n=0 

Generally, one has 

Str(aiai) n = N ^ aiC n ~ l (19) 
i=i 

for some n-dependent coefficients dj. In this section, we will explicitly de- 
termine a , ai and a 2 , whilst in Section 7 we will give the full result for 
Str(ai<yi) n in the case where the generators are in the spin-half represen- 
tation of su{2). 

We will show in the following that 

2 2 
a — 1, ai = -rii n - 1), a 2 = — n(n - l)(n - 2)(7n - 1) (20) 
S 45 

so that 

5tr(a iQ;i )" = N\C n — -n(n -l)C n - 1 + ^-n(n-l)(n-2) (In -l)C n ~ 2 + ... 
\ 3 45 

(21) 

Now note that, if one may write 

Str(a iai ) n = DC n (22) 
for some n-independent differential operator D, then from ()18j) 

oo 

Strfiam) = J2 fnDC n = Df{C) (23) 

n=0 

With the result (J21)) . we see that to the first few orders, 

r, „A 2^ d 2 14^ 2 <9 4 8^ <9 3 \ 
D - N { 1 -3 C ^ + ^ C dC^9 C d^ + ---) (24) 

We will use the results (O, d2H) in Section (4.6). 

3.1 Calculation by evaluation on highest weight 

The basic object of study in the symmetrised trace which we are considering 
is the element ^ 

Cn = J—y^ a ii a h a i2 a h ■ ■ ■ a in a in ( 25 ) 



in the universal enveloping algebra of su(2). The sum is over the (2n)\ permu- 
tations of the a's. Since all the indices are contracted, this element commutes 
with all the generators of su(2), i.e belongs to the centre. By Schur's Lemma 
it is proportional to the identity in any irreducible representation. We can 
calculate it by obtaining its value on the highest weight state of a spin J rep- 
resentation. This method of evaluation on highest weights is quite effective 
in producing formulae for the large J limit of C n , for any n, and is also used 
in Section 7. 

Let us first set down some conventions and useful facts about su(2): 

[cti, atj] = 2ie ijk a k 
a± — a,\ ± 2«2 
[a 3 , a±] = ±2a± 
[a+, a_] = 2a 3 
diOLi = (a^) 2 + + 

The value of a 3 on the highest weight of the spin J representation is 2J 
(where J is 0, |, ■ • • ) - 

a 3 \J, J >= 2J\J, J > (26) 

The quadratic Casimir takes the value C = 4J(J + 1). Evaluation of C n is 
done by writing out the (2n)! permutations and taking the contracted indices 
to be equal to (3, 3) or (+, — ) or (— , +). At the end of this process we have 
a series of a's including a 3 or a±, an example of which is 

(..)a + (..)a_(..)a + (..)a_(..)|J,J> (27) 

The (..) indicate powers of a 3 . We will be more explicit later. All the powers 
of a3 can be commuted to the left say, by using 

a + (a 3 Y = (a 3 - 2) 7 a+ 

a„(a 3 ) 7 = (a 3 + 2) / a+ (28) 

After commuting the a 3 to the left we have an expression made of a string 
of a + and a_ operators acting on the highest weight state. It is useful to 
calculate 

a+c^| J, J >= N(L, m)a L T m \J 1 J > (29) 

where N(L,m) = 2 m ( L ^ m y ^jzxj! • F° r example, for the pattern in (|2*7j) 
we have to evaluate J, J > which is equal to N(l, l) 2 . 
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For L = m = k, the N-factor behaves at large J as N(k,k) = (2J) h . 
Consideration of the iV-factors, together with the explicit powers of 2 J from 
the «3 shows that the highest power contributed by a pattern containing k 
pairs of a± is (2J) 2n ~ k . Hence the C n has an expansion of the form (2 J) 2n (6 + 
bi (2 J) -1 + b 2 (2J)~ 2 + ■ ■ ■ ) and we only need to calculate patterns with at 
most k pairs of a± in order to determine the coefficient bk- 

A useful combinatoric factor needed is the number of arrangements (I25jl 
followed by choices of the values for the i indices, which lead to a fixed (H — ) 
pattern of type (|2~7jl . The factor is fc!( - 2 "~f^ ! | 2 n! . Absorbing the ^yr from 
normalization of the symmetriser we define 

s k\(2n-2k)\2 k n\ . x 

C( *'" ) = (n-t)l(L)l (30) 

The k copies of a + can arise from specifying Oj's contracted to any of k 
copies of Oj's which are specified to a_, hence the k\. The 2n — 2k factors 
of «3 can be contracted with (2n — 2k)\ possible orderings. The first a± 
contraction can be done with one of the n pairs of contracted indices and 
there is a factor of two because we can assign each index of the pair to a 
+ or a — . Hence we have 2n choices for the first contraction, 2(n — 1) for 
the second, 2{n — 2) for the third and so on. Collecting everything we have 
x k\(2n-2k)\ x (2n)(2n-2)(2n-4)... = ^few£^ as claimed above. 
Thus (jHO)) is the factor which must multiply the number obtained by explicit 
evaluation of the (H — ) patterns. 

Let us first consider k = 0. All the terms in the equation defining C n in 
(123)) are equal to (a 3 ) 2n = (2J) 2n . C(0, n) = 1 so the leading term in the 
large 2J expansion is just (2J) 2n 

Now consider k = 1. Let the a_ sit in the ii'th position after c^ 1-1 and 
let the a + sit in the ^'th position after an additional ag 2- ' 1-1 



l 3 

2n-l 2n 



a 3 a+ctg a_a 3 \J,J> (31) 

h=l i2=ii+l 



Multiplying by the combinatoric factor in (p?U|) and commuting the a.3 to the 
left we get 



2n-l 2n 



C(l,n)J2 "3 n "' 2 («3 -2) i2 ~ il - 1 a^ 1 a+a_|J, J> 

h=l 12=11+1 
2n-l 2n 

= C(l,n)J2 (2J) 2n ^ 2 (2^-2) i2 " n - 1 (2J) il - 1 iV(l,l)|J, J> 



il=l 12=H+1 
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which can be simplified to 



= C(l,n)N(l, 1)(2J) 2 ^ 2 ^(1 - \r- h - 

i 

= (2J) 2 - 1 ( 2n - \n{n - 1)-L + \n{n - l)(2n - 3)- 



2J 3 v /v 7 (2J) 2 
1 fi 1 
-_n(n-l)(n-2)(2n-3)^3 ) (32) 

The sums here and below can be done with a mathematical software package 
such as Maple. For k = 2, there are two patterns : 

(..)«+(••)«-(••)«+(••)«-(••) (33) 
(..)«+(..)«+(••)«-(••)«-(••) (34) 

We will denote them as the (H 1 — ) pattern and the (+ H ) pattern. 

An alternative notation to distinguish them is to write 

1 1 
1 1 

for the first, and 

'2 

for the second. In the first array, the first integer in the first column indicates 
the number of successive a_ seen while we read from the right and the second 
gives the number of a + that follows. The top integer in the second column 
gives the number of a_ that follows after this and the lower the number a + 
thereafter. In the second array, the upper 2 is the number of a_ seen reading 
from the right, and the lower 2 is the number of a + after that. Steps similar 
to the case k — 1 give for the value of (+ H ) on the highest weight 

2n-3 2n-2 2n-l 2n 

7V(2,2)C(2,n)(2J) 2 "- 4 J] E E E (1 - -) i2+ — i3 ' 2 (35) 

11=1 «2=«1 + 1 «3=*2+l «4=i3 + l 

Evaluating the pattern (H 1 — ) gives 

iV(l,l) 2 C(2,n)(2J) 2n - 4 

2n-3 2n-2 2n-l 2n 

x E E E E (! - 7 )- +i2 - ii ~ i3 - 2 (i - 7 ) 1 — - 1 (36) 

11 = 1 12=U+1 i3=*2 + l 44=13+1 
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These sums can be evaluated and collecting the contributions for k = 2 we 

get 

1 



(2J) 2n " 2 ( 4n{n - 1) - -n(n - l)(4n - 7)- 



3 v /v 7 (2J) 

X 6 X 

+ - l)( n - 2)(50n - 101)^ ) ( 37 ) 

For = 3 we have 5 patterns which we list below together with the 
corresponding iV-factors - 



3 
3 

2 1 

2 1 

2 1 

1 2 

1 2 

1 2 

1 1 1 

1 1 1 



7V(3,3) 

7V(2,2)7V(1,1) 
/V(2,l)/V(2,2) 
/V(l,l)/V(2,2) 
JV(1,1) 3 



(38) 



Using the compact notation J2i f° r 

2n-5 2n-4 2n-3 2n-2 2n-l 2n 

££££££ 

»1=1 12=H+1 i3=«2 + l «4=*3 + l «5=*4+l «6=*5 + l 

the result of evaluating the five patterns above is 

N(3, 3)C(3, n)(2J) 2ri ~ 6 ^(1 — — ) i2+i «- il - i5 - 2 (l— - ^-Hs-ia-u^. - y^-h-i 

(40) 

iV(2,2)A^(l,l)C(3,n)(2J) 2ri - 6 ^(l - I^+^+w-w-w-is-s^ _ )is-«-i 

(41) 

N(2,l)iV(2,2)<7(3,n)(2J) 2 "- 6 ^(l--)^ 

(42) 
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iV(l,l)iV(2,2)C(3,n)(2J) 2n - 6 ^(l - I^+^-n-^-s^ _ 2y B -u-i 

'J U 

(43) 

N(l, 1) 3 C(3, n)(2J) 2n - Q ~ -i)* 2 ^-^ 6 -* 1 -* 3 -* 8 - 3 (44) 

i 

After expanding in 1/ J and doing the sums, and collecting the five con- 
tributions for k = 3 we get 

(2J) 2 "- 3 (8n(n - l)( n - 2) - 16(n)(n - l)(n - 2)(2n - 5)-^) (45) 

For k = A there are 14 patterns. We write the array description of the 
patterns, followed by the corresponding iV-factors - 



4 

3 1 

1 3 

3 1 

3 1 

2 2 
2 2 

2 1 1 

2 1 1 

2 2 

1 3 

2 1 1 

1 2 1 



iV(4,4) 

7V(3,3)7V(3,1) 

7V(3,3)7V(1,1) 

iV(2,2) 2 

iV(l,l) 2 iV(2,2) 

7V(3,3)7V(2,1) 

iV(2,l)iV(2,2)iV(l,l) (46) 
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and 



2 1 1 
112 

1 3 
1 3 

1111 
1111 

1 2 1 
112 

112 
112 

1 2 1 
1 2 1 



N(2, l) 2 iV(2,2) 

iV(l,l)JV(3,3) 

AT(1,1) 4 

jV(l,l)JV(2,l)iV(2,2) 
N(l, l) 2 iV(2,2) 

iV(l, l) 2 iV(2,2) (47) 



For each of the 14 terms there is a sum of the type in ()39|) generalised to 
8 summation variables. After (2J) 2ri ~ 8 C(4, n) multiplied by the appropriate 
iV-factors is extracted, what is left is nothing but the number of terms in 
the sum which is ^ 2 n-sy.8\ • Collecting all patterns relevant to n = 4, with the 
appropriate combinatoric factors, and evaluating on the highest weight gives 

16(n)(n - l)(n - 2){n - 3)(2J) 2 ™- 4 (48) 

With the results in equations (|32|37|45|48J) . we can find the value of C n 
SIS db function of J to the first few orders 

C n = (2J) 2n (l + (2n)^- + U{n-l)-^-,-U{n-l) ' 



2 J 3 v J {2J) 2 3 v J {2J) 



+ ^n(„_l)(„_2)(n + 2)^_ + ...) (49) 

Matching this with an expression of the form aoC n + aiC n ~ l + g^C" 1-2 + • • • 
where C = 4 J( J + 1) determines 

a = 1 

2 

ai = --n(n-l) 

a 2 = — n(n-l)(n-2)(7n- 1) (50) 

It is worth noting that once we have determined the coefficient a® by 
calculating the patterns with k = which determine the coefficient of (2 J) 2n , 
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the formula in terms of Casimirs fixes the order (2 J) 2n_1 term. The latter is 
checked by considering patterns with k — 1. Further consideration of k = 2 
fixes the coefficient a\. With ao and a x fixed the term of order (2J) 2n ~ 3 
is fixed by the expansion in powers of the Casimir and can be checked by 
considering patterns with k = 3. Expanding the contribution from k = 3 
patterns to next order and collecting the leading order from k = 4 allows the 
determination of ci2- With ai,a2 determined, we now have a prediction for 
the coefficient of (2J) 2n ~ 5 , i.e ^n(n — l)(n — 2){n + 2). We have checked by 
considering the 42 patterns which arise at k = 5 that this is indeed the correct 
coefficient. Independent confirmation of these results using computations 
based on evaluation of the Casimirs using chord diagrams are given in the 
following. 

3.2 Casimirs and chords 

The group invariants discussed above arise in knot theory, in particular in the 
study of finite type invariants of knots. These form a type of basis of knot 
invariants, and may be understood as terms from the perturbative expan- 
sion of Chern-Simons gauge theoretic knot invariants. Finite type invariants 
may be discussed in terms of chord diagrams [21], which are diagrammatic 
representations of precisely the individual terms which occur in Str(a;aj) n . 
A term in this symmetrised trace will be the trace of a product containing 
n pairs of group generators in some order a^a^-.-.a^. Writing each 
matrix generator as a vertex 



the matrix product joins matrix-labelled legs of the vertices, with the trace 
forming a circle. Chords of the circle are then formed by the pairs of free 
legs with like indices. For example, 



(conventionally one starts at the point on the circle corresponding to twelve 
o'clock and moves counter-clockwise around the circle as one moves from left 
to right inside the trace). An individual term in Strict;)™ will be represented 
by a circle containing n chords, ie a chord diagram. Particular choices of 
groups and representations then give explicit realizations of chord diagrams, 




a 



b 




14 



the resulting assignment of numerical values to diagrams leading to a "weight 
system" and a finite type invariant of knots. In the following, the "order" 
of a chord diagram will be the number of chords, and the "value" < D > of 
a chord diagram D will be the numerical value obtained by evaluating the 
group theoretic trace for the group and representation under consideration 
(for simplicity, we will ignore the factor coming from the trace of the identity 
matrix in the expressions for chord diagrams in this section). 

An operator of much interest for knot theoretic applications is that of 
cabling. This is closely connected with the Adams operation in group theory 
and with the fundamental Alexander- Conway invariant in knot theory [22] . 
Cabling involves taking multiple covers of the encompassing circle in chord 
diagrams, and lifting the chord ends in all possible ways to this new circle. It 
is clear by definition that the symmetrised trace Str(ajO!j) n is invariant (up to 
a factor) under the cabling operation. In fact it is the unique eigenvector of 
highest eigenvalue of the transpose of the cabling matrix [22] • For example, 
for the chord diagrams with two chords, define the basis 




Then the cabling operation ■0, which takes the double cover of the chord 
diagram circle and lifts chord ends, is given by 

* = *(a a) <«) 

The eigenvectors of ip T are then (—1, 1), (2, 1), with eigenvalues 4, 16 respec- 
tively. For comparison, Str^a,) 2 = |(2, 1), which is proportional to the 
eigenvector with highest eigenvalue. We note that the symmetrised trace 
STr includes a normalization factor l/(2n)!. 

There are in general relations between chord diagrams [21], arising from 
identities satisfied by the generators aj. For example, at order three, 




Taking the full set of such relations into account, one can show that the 
following sets of elements define bases for diagrams of order 2, 3 and 4 - 
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We will now specialise to the group sl(2). In this case, it is possible to derive 
a reductive formula relating the value of a chord diagram with n chords to 
the values of combinations of chord diagrams with fewer than n chords 
(the group theory conventions of this reference will be used in this section. 
This will introduce some factors of 2 as will be noted below). Thus one 
may express the value of any chord diagram as a polynomial in c, where the 
Casimir c is the value of the diagram with one chord: 




(52) 



Using this reduction to evaluate the diagrams in the basis elements defined 
above, one finds the explicit results 

(c 2 , c(c - 2)) 

(c 3 ,c(c 2 -6c + 8),c(c-2) 2 ) 

(c 2 (c - 2) 2 , c(c - 2) 3 , c(c 3 - 10c 2 + 40c - 40), c(c - 2) 2 (c - 4), c 3 (c - 2), c 4 ) 

The expressions for Str(aiaj) n , for n up to 4, may then be found, either by di- 
rectly carrying out the symmetrised trace prescription, or by using the action 
of the cabling operator ip and deducing the highest eigenvalue eigenvector of 



16 



i\) T . One finds for general groups 

Str(ajCtj) 2 
Str(a i aj) 3 
Str(ajaj) 4 

Evaluating these for si (2) yields 

Str^tt;) 2 | s z(2) 

Str(a i a i ) 3 | si ( 2 ) 
Str(ttjaj) 4 | si (2) 

Evaluating a chord diagram gives a polynomial in c. The coefficients of this 
polynomial are given in terms of quantities defined by chord intersection 
properties of the diagram [22] • For a chord diagram D n with n chords, one 
finds 

< D n >= c n - 2Ic n - 1 + (21(1 - 1) - AT + 8Q)c"- 2 + o(c"- 3 ) (53) 

where I is the number of intersections of pairs of chords in the diagram, T 
the number of triple intersections and Q the number of quartic intersections 
(intersections of four chords in the shape of a square, with possible other 
intersections amongst these four chords). This leads to corresponding results 
for the symmetrised trace operator. Denote by < X > av (n) the value of some 
quantity X defined for chord diagrams, averaged over the set of all chord 
diagrams of order n which arise from the symmetrised trace. Then, for the 
result quoted above we see that the coefficient of c n_1 in the polynomial 
arising from evaluating Str(aiai) n is just —2<I > a v(n)- 
We have 

< / >av(n)= ^n(n - 1) (54) 

To prove this, note that any pair of chords can either intersect or not in any 
diagram. Thus, to find < I > ov ( n ), we consider all possible positions of a pair 
of chords on a circle where there are 2n possible endpoints for chords. For 
each such placement of the pair, one sees if / is or 1, and then sums this, 



|(2,D 
^(1,-1,3) 
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(6,2,1,1,3,2) 



c c 

3 



c A - 2c 2 + -c 



c 4 - 4c 3 + — c 2 - — ( 
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and finally one divides by the total number of such placements. This leads 
to the expression for each pair of chords 



Y^\(2n-2-i) _1 



(n- l)(2n- l)(2n- 3) 3 

This is to be multiplied by the number of choices of pairs of chords, ie |n(n — 
1), giving the result above. Alternatively, since the result is proportional 
to n{n — 1), one can find the coefficient by computation for n = 2. This 
provides a different proof of the result in Section (2.1) above (to compare 
coefficients, we note that the Killing form used in j2H| introduces a factor 
of 2~ % multiplying the coefficient of c n_? , and so one needs to divide by this 
factor in translating to the conventions used elsewhere in this paper). 

For the next order term, involving c n ~ 2 , to calculate < T > a v(n) and 
< Q >av(n), note that these are proportional to n{n — l)(n — 2) and n(n — 
l)(n — 2)(n — 3) respectively. The coefficients are fixed by explicit calculation 
for n = 3, 4 respectively. This gives the results 

< T > av{n) = — n(n- l)(n-2) 

yu 

< Q >av(n) = Y^n(n - l){n - 2){n - 3) 

For < I 2 >av(n): note that at large n this behaves as n 4 , and it vanishes for 
n = 0,1. This fixes all but three of the coefficients of the expression as a 
fourth order polynomial in n. These further three coefficients are then found 
by calculation for n — 2, 3, 4. Thus one finds that 

</ 2 >av(n)= ^-n{n - \){hn 2 -n + 12) (56) 
loU 

Putting together the above expressions, inserting the factor of iV coming 
from the trace of the identity, and using the conventions of the rest of this 
paper, we are thus led to the result 



Str(a iai ) n U(2) = C n --n(n-l)C n - 1 +—n(n-l)(n-2)(ln-l)C n - 2 +o{C n ~^ 



1 

_ ;v ' '45 

(57) 

in agreement with the calculations of the previous section. 

Clearly the precise form of subsequent terms in this expansion will depend 
upon detailed features of diagrams of increasing complexity. In order to 
obtain information on the general structure of these terms, we now outline 
an argument concerning the large n behaviour of Str(ajaj) n . First note that 
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at large n, ie for diagrams with a large number of chords, the symmetrised 
trace generates sets of diagrams which are dominated by those with / large 
and the other quantities T, Q, etc, relatively small, since / is always greater 
than these quantities in any given diagram, and there are in addition many 
diagrams with I non-zero and the other quantities all zero. To determine 
the pure / dependence of a chord diagram, consider the diagram P n which 
corresponds to the expression 

Tr ^(aj 1 aj 2 ai 3 aj 1 )(o;j 4 Q;j 3 aj 5 aj 4 ) ( a «2n-3 Q; *2n-4 a «2n-2 a «2n-3) 

X ( a i2n-3 a hn-i a i2n-2 a i2n a i2n-i a i2n a i2)^J (^8) 

This corresponds to a chord diagram, with n chords in the shape of an n- 
polygon inside the circle, but with one pair of adjacent chord ends (labelled 
by a 2 and a^ n in this case) interchanged to remove one intersection. The 
diagram P n has I non-zero and all other geometric quantities T, Q, etc, zero, 
and so it can be used to deduce the pure / dependence of a general diagram. 
It is not difficult to prove by calculation that the value of P n is given by 

< p n >= c (c - 2)™- 1 (59) 

This indicates that the pure / dependence of a general diagram D n with n 
chords is given by 

< D n >= c n - 2ic n ~ 1 + 21(1 - l)c n ~ 2 + • • • + (-2) l c n - 1 + ... 

+ ( n -i)<- 2 >"" c < 6 °) 

We now wish to average this over chord diagrams at order n. Following 
arguments similar to those given earlier, we find that at large n we have 

< I x > av(n) ~ (n 2 /6) x (61) 

Using (joTIj) . we can deduce the contributions to the symmetrised trace from 
chord diagrams with up to double intersections at large n. This gives an 
approximation to the symmetrised trace which is (in the conventions of the 
rest of this paper) 

-str^r - Ei! (t) n2l ° n ~ l 

i=o ' \ ' 

= C n - -n 2 ^- 1 + -n A C n ~ 2 - — n 6 C n ~ 3 + . . . (62) 
3 9 81 V ' 
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The expression (|fi2*|) shows that the signs of terms of decreasing order in 
C alternate, and that the coefficients fall to zero somewhat faster than an 
inverse factorial. This formula gives a guide to the general structure of higher 
order terms. From (JSTj) the large n limit of the exact expression is C n — 
(2/3)n 2 C n ~ 1 + (14/45)n 4 C n ~ 2 + . . . , and we see that ignoring the effects of 
higher order intersections in diagrams has reduced the coefficient of the C n ~ 2 
term by about 25%. The effects of higher order corrections with the structure 
in (pT2*|) will be reviewed in Section 6. 



4 Energy corrections to first order in 1/N 

In this section we will study the effects of the corrections to the theory defined 
by (fT5^) which are of the lowest order in 1/N. The leading order Lagrangian 
is 



Co = -Vl +4A 2 C7? 4 V1 - X 2 CR 2 (63) 
and the energy to lowest order is 

En = R : £n 

OR 



'l + AX 2 CR 4 
1 - X 2 CR 2 



|— (64) 
1 — s z 

In the last line we used the dimensionless variables in ()13j) . Including the 
corrections to next order, by the arguments of Section 3 we are led to the 
corrected Lagrangian 

/ 2 d 2 \ 

c > = I 1 - -f ocp) c " (65) 

We will work with re-scaled variables r 4 = 4X 2 CR A and s 2 = X 2 CR 2 . The 
conserved energy to this order is then 

£ ^*jg- c ^{ l -l c &) £ ° (66) 

It is useful to use the variables U and 7, defined by 

U = Vl + 7\ 7 = ,} 2 (67) 



VI - 
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Then the zeroth order energy is just 

£o = jU (68) 

which indicates that one might think of U as a position dependent "mass" . 

We will denote the first order expression for the energy by the symbol E 
in this section. From equation ()66|) this is given by 

E(j, U) =1 U- ^(l 2 U 2 - l)(3 7 2 f/ 2 - + 1) (69) 
4.1 Plots 

It is instructive to study the energy (p)9"|) as a function of the variables (r, s), 
using (jHZj). We will do this in the following. 




r 



Figure 1: Plot of the energy as a function of r, for s = and taking C = 50. 

Plots of the energy function as a function of r show that at s = it 
is monotonically increasing. Figure ^ shows such a plot. This behaviour 
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means that as the D2-brane starts collapsing from rest at some large radius 
ro, there is no smaller radius where the energy can take the same value at 
t = while returning to zero velocity. If the potential had a minimum, 
a conventional bounce would be possible with the radial evolution slowing 
down to zero velocity and then re-expansion starting with a change in sign 
of the radial velocity. So clearly such a bounce does not happen with the 
first 1/N correction. 
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Figure 2: Plot of energy as a function of s, for r = and with C = 50. 

Plots of the energy as a function of s at fixed r show an extremum. Figure 
|2] shows such a plot at r = 0. If the energy is bounded above at r = 0, this 
means that a D2-brane, starting from zero velocity at sufficiently large r$ 
(hence sufficiently large energy), cannot reach zero radius. This may seem 
paradoxical given the conclusions to be drawn from the plot above of energy 
as a function of r at s = 0. Further insight can however be obtained by 
plotting contours of constant energy in the (r, s) plane. 
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Figure 3: Plot of constant energy contours as a function of r and s, with 
C = 100. The energy of the contours is increasing from left to right. 
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Figure El shows contour plots of the energy (fHHJ) as a function of r and 
s. We have chosen C = 100 as an example. If we do not include the \ jC 
correction the constant energy contours all reach down to zero radius. The 
key feature is that ^ = at a point where s = s m ^ 0. This does not happen 
for ordinary mechanical systems described by energy functions of the form 
E = \ms 2 + V(r). For such energy functions fjf = only at the point s = 0. 
As we will see from analytic considerations in section 4.2, solutions to ^ = 
or equivalently |p = do exist for 7 ~ C 1//4 . 

More detailed numerical investigation of the trajectories are possible. For 
fixed initial energy Eq or initial position tq (or equivalently Uq) where the 
velocity is zero, we can solve the equation E(j,U) = E for U(j,E ). The 
equation is quartic and has four solutions, one of which is the physical one. 
This solution can be followed as 7 is changed from unity to some large num- 
ber. For Eq corresponding to Uq > 1 but Uq <~ C 1 / 4 we find that the 
classical path describes U decreasing monotonically to U = 1. These paths, 
which we may call perturbative paths , are small ^ perturbations of the paths 
obtained from Cq. For larger values of Uq, the path encounters a minimum 
radius at some finite value 7 m and proceeds to increasing U, which in fact 
approaches infinity at a finite upper bound. 

Figure 0] shows the classical path in the (U, 7) plane for C = 50 4 and 
Eq = 20. For these choices, U decreases continuously through U = 1 or 
r = 0. Figure El shows the classical path in the (U, 7) plane for C = 50 4 
and Eq = 100 2 . For these choices, U decreases until it reaches a minimum 
and then increases apparently to infinity at a finite value of 7. Beyond this 
limiting 7 the quartic equation for U given by E(j, U) = Eq has no physical 
solution. Along these paths we can also compute an effective mass \J E 2 — p 2 
and the (proper) acceleration. We will discuss these further in Section 8. 



4.2 Extrema 

The extrema in the energy contour are one of the novel features of C\. As 
we will elaborate on later, they occur in a region where the 1/N expansion 
is becoming problematic, but we will nevertheless study them in more detail 
since it is entirely possible that such extrema occur in the finite N case. 

On a contour of constant energy one has dE = 0. To first order in 
variations of 7, U we can write dE = (^p)u^7 + (|^) 7 dC7. At an extremum 
of position along the energy contour, we have dU = 0. This means that 

(i% = (70) 
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Figure 5: Plot of U as a function of 7 = g along a solution, with C = 50 and 
E = 100 2 
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at that point. Equivalently we have 



®7 1 JTT ' - 



)E 



(71) 



If we consider U as a function of 7 along a constant energy contour, an 
extremum of U requires the vanishing of the energy derivative (|70j). 
This gives us a quadratic equation for U 2 in terms of 7 

This is solved by 



= 2(-l+5^-¥) 

For 7 close to 1 the denominator is negative, while the argument of the 
square root is negative as it is dominated by ^ at large C, so there is no 
real solution for U. The positive branch must be chosen for a real solution 
to be possible. For C large compared to 1, and 7 2 ~ y/C, the denominator 
changes sign and there is a real solution. The minimum 7 which allows an 
extremum is obtained by setting the denominator to zero. This gives 

7 2 = l -{2 + V4 + IOC) ~ ^VC (74) 

where we have used the large C approximation in the last step. The associ- 
ated minimal velocity is s 2 = 1 — 

When 7 is very large compared to C, we again have no physical solution, 
because U becomes close to zero, whereas for real r it must be larger than 1. 
The upper bound on 7 for allowed extrema is obtained by setting U — 1 in 
(I72|) and solving for 7. This gives 



7 2=I( 3 + V4 + 10C)~y^ (75) 

Comparing the above two equations, we see that there is a very thin range 
of 7 which allows the extremum to exist. 

While there is a finite range of velocities which allows the extremum to 
exist, there is no bound on r. This is a bit surprising since we expect the 
extremum to occur in a sense near the origin. This can still be true in the 
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following sense. If we start the collapse with zero velocity at some position 
Uq, and ask for the minimum U reached, we will get a minimum which is a 
function of U . Let us denote this minimum as U m (U ). We would like to 
prove that ~ ± where a is a positive number, so that this approaches 

zero in the large C limit. Now consider a brane starting with zero velocity. 
Substitute 7 = 1 in (ffiTJj) to obtain the energy E at the starting position Uq 

E (Uo) = U Q + -^(ut-l) 2 (76) 

We can determine 7 as a function of U along the curve of extrema by solving 
the quadratic equation for 7 2 in f72J) . This gives 

2 = (3 + 6gy + y/24 - 24Uj + gg + 90£%C 

Substitute this in the energy formula (jBTJj) to get a function E(U m ) given by 
2(45 + 90££ + 15D)3(90C^C + 18 - 48£/£ + D + 12E# + 1U 2 m D) 



E(U n 



(78) 



3375C/ 4 C 
where D is given by 

D = V(24 - 24£ft + 36tft + 90tftC) (79) 



Equating E{U m ) = E(Uq) gives an equation which determines U m (U ). This 
becomes tractable if we use the large C approximation (assuming U m is 
sufficiently small compared to C). We find that U m /Uo = C -1 ' 4 /o(^o) + 
C- 3/4 /i(f/o) + • ■ ■ , where f (U ) ~ 1 + 0(£) and /^C/q) ~ 1 + 0(£) ■ 

4.3 Energy-momentum relation 

We now turn to dispersion relations. The momentum is 

1 dL 



P 



\VC ds 



= 7^-^(3t/V + l)(^V-l) (80) 

We have included a factor of —7=? for convenience. Note that the formula 
for the momentum, like that for the energy, has at the leading term the 
standard form of special relativity, with U playing the role of mass. We 
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would like to obtain an equation for E(jp) which can be interpreted as a 
deformed dispersion relation. 
Let the corrected speed be 

P (81) 



vV + u 2 c 

where we will determine s\ as a function of U and p. The corrected 7 variable 



is 

7 



Substituting this in ([5UJ) and keeping the leading term in the 1/C expansion 
gives us a linear equation for s±. The solution of this is 

3i = lP ^ J (3f/ 2 + 3p 2 + l)(f/ 2 + P 2 - l)p (83) 

Substituting in the corrected speed (or the corresponding corrected 7) in (|69j) 
we find the first-order corrected energy formula as a function of p is 

E = yfp^IP + ip2+ J^ 3/ \ u 2 +p 2 -l)(U*-3p*-l)(U 2 + P y 
2 

P (3U 2 + 3p 2 + 1)(U 2 +p 2 - 1) (84) 



6CU 4 

There has been a lot of discussion of deformed dispersion relations recently 
in the context of inflation (see [21], for example). The deformed dispersion 
relation we are getting does not have a smooth limit in the zero "mass" limit 
of U = 0. In our context this limit is of course unphysical since U = a/1 + r 4 . 

4.4 The time of collapse 

It is also possible to determine the time of collapse in this system. The 
leading formulae for velocity s and the corresponding 7 are 



S(o)(U, U ) = U \ 7(o)(t/, U ) = ^ (85) 

We will find the perturbed formulae 
* = *(o) + ^s {1) (U, U ), 7 = 7(0) + ^7(o)(l - 7(o))* = 7(0) + lj §- (86) 



29 



where the last line defines 7(1). We substitute fjTfij) into the left-hand side of 
(jnilj), and on the right-hand side put in the 1/C expansion of the velocity, 
keeping only the leading terms. This gives 

E (U ) -U = ^ + i/( 7( o), U) (87) 

where / is the correction term that appears in (jfffljh This allows us to solve 
for 7(x), or equivalently S(i). We find 



HD = ^iu^^^ ^V o - UlU 2 + Ul + f/ 2 ) (88) 



Thus the 1/C corrected formula for the speed is 



A" ^ + C~6^ ^ (3U +U +U (l-U )) (89) 

and the time of collapse is 



These integrals can be evaluated in terms of elliptic functions. 



5 Solution in the neighborhood of the ex- 
tremum 

In this section, we will study the behaviour of the brane in the region of 
the radial extremum found above. Generally, consider the expansion of the 
energy around the extremum, keeping the second order terms. Consider 
a constant energy contour parameterised by A, i.e we think of coordinates 
(G(X),V(X)), where G is some function of the f, for example 7 or f. V 
is some function of r, for example U or more simply r. By expanding the 
energy in a Taylor series we find at lowest two orders 

„fdEdV dEdG\ n 

HwiX + dGlx)=° (91) 
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and 

J\dEd 2 V ldEd 2 G l/dVyd 2 E 1 /dGyd 2 E d 2 EdGdV 
\2Wli^ + 2dG~d)^ + 2\~d\) dV 2 + 2\~d\) dG 2 + dGdV~d\~d\ 

(92) 

Specializing to a point where = 0, then from (j9~T|) we have §§ = 0. Then 
we obtain 

ldEd 2 V l/dVyd 2 E l/dGyd 2 E d 2 E dGdV _ 
2dV~d>? + 2\d\) W 2 + 2\1\) dG 2 + dGdVlX ~d\ ~ ° ^ 



If we approximate the equation in the neighborhood of the extremum where 

dV 
<IX 



4j = by setting this first derivative to zero we get 



_ dG^dx) (qA) 

The complete equation derived at order d\ 2 is really (|9l?|) but we may hope 
that some qualitative properties such as the existence of two branches cor- 
responding to reflection with increasing speed or reflection with decreasing 
speed can be captured by the simpler equation (jHlj). This equation should 
allow us to prove that the extremum is a minimum of U and that we need a 
change of branch of the solution at the extremum. 

Let us choose V = r,G = s 2 = S and A = S so that (JiMj) becomes 

d 2 r , , 

w = -a. P5) 

where a = |§^f7fjf evaluated at the extremum. Equation (J93|) can be solved 
in the neighborhood of the minimum to give 

(r-r m ) = -a(s 2 -s 2 J 2 (96) 

We know that r > r m because plots show that the constant energy contours 
have a minimum and a is less than zero. We can now write the squared 
velocity 

s 2 = s 2 m ± sj^- (97) 

Hence 



IrWW 1 ^ (98) 
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which can be solved for the time elapsed 



/ 



dt = ± / dr 




1 



(99) 



a 



When we choose the outer negative sign, r decreases with increasing time, 
whereas when we choose the outer positive sign, r increases with increas- 
ing time. The collapsing solution cannot be continued past r m because the 
formula for time elapsed would become complex. These two solutions can 
be patched at r = r m . The radial velocity is discontinuous at the patching, 
but the time evolution described by the patched solutions is consistent with 
energy conservation. 

Note that there is also a choice of sign inside the square root. If we 
switch the choice of this sign in patching solutions at the extremum, then 
the collapsing brane re-expands along the second branch, which means that, 
when it reaches the original starting position (where it had zero velocity), it 
is travelling close to the speed of light. If we do not make this switch of sign 
inside the square root upon patching, the brane re-expands along the same 
branch as it was collapsing and reaches zero velocity at the starting position. 

If we use the form (V, G) = (s, r) and use A = s, the differential equation 
flUH becomes 



Here the overall choice of sign visible in (}9~9"j) is not apparent. However we 
know by the time reversal invariance of the Lagrangian that solutions which 
have an extra ± in front of the integral are also allowed. This time reversal 
invariance was kept manifest when we worked with the variables (r, s 2 ). 

A more precise treatment of the local differential equation in the neigh- 
bourhood of the extremum continues to reveal that s — s m is double valued 
there. With the choice (G, V) = (s,r), A = s, if we keep terms involving the 
first derivative we get an equation of the form 




(100) 



where f3 = |0/§f . Then 




(101) 




(102) 
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where the constants are given by 



, fo =^-' ^=^T ( 103 ) 

Equation (jl()2j) is the Riccati equation with constant coefficients. We can 
solve this as follows. Use variables y = j- to rewrite (|102j) as 

a^ + by 2 + cy + d = (104) 
as 

The roots of the quadratic polynomial in y play an important role in the 
solutions. These roots are 

-c + D -c-D . , 

« + = a_ = (105) 



where we defined .D = \J c 2 — Adb. Numerical check shows that D is real. 
Integrating once gives 

ds=-l r J—An(^±) (106) 
b [a + — a_) \y — a_ / 

Imposing the condition that y = ^ = 0ats = s m fixes the integration 
constant to give 

<— > = £ (-e)->"0) 

One more integration and the condition that r = r m at s = s m gives 

. a ( a + fb(a_ - a + ) 
r - r m = a + (s - s m ) + T In 1 exp (s - s„ 



6 \ a_ V a 

- H 1 -^) (io8) 

The original differential equation (|104|) is symmetric under exchange of the 
roots a± and as expected the final solution can be manipulated into the form 

. , a ( a_ ( b(a + — a) 
r - r m = a_(s - s m ) + —mil exp (s - s 



b \ a + \ a 
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Expanding the right-hand side of (jl()8j) or gives 



(r - r m ) 




) 



2 



(110) 



2E i r m) 



s 



s 



m 



in agreement with the approximation which dropped the j- terms in the 
equation. From (|11U|) it is clear that for a fixed r there are two values of s, 
one above s m and one below. 

We can also use the time t as the parameter A and obtain an equation 
which looks a little more complicated but can nevertheless be solved. It seems 
that the time variable may be less useful since dr/dt has a discontinuity at 
the extremum. In any case the equation is now 



This can again be solved. 

5.1 Quantum mechanics near the extremum 

Equation (jl(J(J|) can be viewed as describing the evolution of r as a function 
of a "time" variable s. We can write down a Lagrangian which leads to 
the equation and a corresponding Hamiltonian which generates translations 
along the s direction 



The solutions to this equation are Airy functions. The correct boundary 
condition is to require that the solutions die off for r < r m . This implies 
oscillatory behaviour with a damped amplitude as the magnitude of s — s m 
increases. It is important to note that in this quantum mechanical set-up 
s < s m and s > s m are treated symmetrically. 

When the terms are kept the system becomes dissipative. Such systems 
cannot be quantised in an ordinary way but they can be described quantum 
mechanically by enlarging the system to include a bath of bosonic oscillators, 
see for example |2*5| EH] . ( A connection between string theory and this 
quantization of dissipative systems has in fact been made in [27] ). The 
correlation functions < r(s)r(s' > can be considered. Such treatments in 
terms of a bath of oscillators can be viewed as simulating the interaction of 
the massless degrees of freedom including r with the stringy spectrum of open 
string modes. It is possible that the extra degrees of freedom required are 
just the additional modes that live on the brane world-volume, e.g the various 





d 2 i) 



+ brip = Eifj 



(112) 



dr 2 
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fluctuations of the fields (which have been neglected in the solutions where we 
study purely radial evolution). The important conclusion for our purposes 
is that whereas the classical evolution is ambiguous at the turning point, 
a quantum mechanical set-up allows computations of correlators associated 
with different points along the trajectory on either side of the extremum. 
This leads us to anticipate that a more complete stringy treatment will allow 
the formulation of probabilities for the evolution along the two branches. 

It may appear that the conclusions from attempting the above quantum 
mechanical discussion are artefacts of the exotic choice of s as time. However, 
the ordinary time can also be used to describe the equations, which are 
then of higher order. This again suggests that we need additional degrees 
of freedom. In fact it has been argued that particle systems with higher 
derivative actions share some stringy properties such as the Virasoro algebra 
28J. The conclusion from this discussion is again that an attempt to do 
quantum mechanics leads to the need for extra degrees of freedom, which 
are indeed expected in this stringy context. While these stringy modes are 
not relevant for most of the trajectory, at the extremum the only possible 
classical evolutions involve discontinuities in velocity, hence formally infinite 
accelerations. Since small accelerations are a requirement for the low energy 
effective action to be useful, in this case we expect that extra degrees of 
freedom become relevant. 



6 Higher order corrections to the Lagrangian 

Based on the results of Section 3, the next order corrections to the Lagrangian 
of (JfiHj) give the result 

{ l ~ l c & + + l c &) c ° (113) 

The energy function to this order is 

7 f_. TT , i\/..rr i\/o..2rr2 , i att2\ 



£2 = iU--^( 1 U + l)( 1 U-l)(3 1 2 U 2 + l-AU : 



7 



r( 7 C/ - l)(jU + 1) (245 7 6 [/ 6 - 640 7 4 f/ 6 + 105 7 4 £/ 4 



120C 2 U 7 

+ 528 7 2 f/ 6 - 256 7 2 f/ 4 + 63 7 2 £/ 2 - 128f/ 6 + 176£/ 4 - 128£/ 2 + 35 

It is possible to repeat the analysis of Section 4, studying the effects of the 
above terms of order 1/C 2 . We omit the details. In summary, as in the study 
of the Lagrangian £ 1; there is a class of perturbative paths in (7, U) space 
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(for Uq <~ C 1 / 4 ), which are similar to the ones arising from the Lagrangian 
Co- For higher Uo the paths still end at U — 1 without going through 
the extremum seen with C%, but there are important differences. The proper 
acceleration does not increase far beyond 1 and the final 7 at U — 1 increases 
with Uq much slower than would be anticipated by extrapolation from the 
final 7 for the perturbative paths. This suggests some sort of repulsive force. 
It also indicates that there may be a class of non-perturbative (in the sense 
of the 1/N expansion) paths having small proper accelerations which allow 
simple reliable treatment with the low energy effective actions neglecting the 
higher derivative terms. 

Going to one further order, using the term suggested by the expression 
(JH2J), however restores the extremum when these terms contribute. Thus, as 
this emphasises, in regimes where higher order terms are not negligible it is 
necessary to obtain information about the exact structure of Str(ceiai) n . To 
this end, in Section 7 we will derive and study the exact expression for this 
when the spin half representation is used. 

7 Exact evaluation of the symmetrised trace 
for spin half 

In this section, the C n operator will be evaluated exactly for the spin 1/2 
case. We will be using some of the notation and techniques from Section 
(3.1). Consider patterns of the form 

al^a+ai 2 "'^- 1 ' 1 ■ ■ ■ a^ 2 " V^ 2 " 11 " 1 "-^ 1 " 1 (114) 

These are the only patterns that contribute for spin 1/2 since any pattern con- 
taining (. .)«_(. .)«_(..) (where (..) as in Section 2 stands for arbitrary powers 
of CK3) gives zero when acting on the highest weight of the two-dimensional 
representation. This is an important simplification compared to the case of 
general J. Expressions of the form ()114j) are to be summed with summation 
symbol J2i defined as 

2n «2fe-l 13—1*2—1 

E E -EE (us) 

22fc=2fc i 2 fc_i=2fc-l 12=2 ii=l 

After commuting all the 0:3 's to the left we get 
VVai3 — 2)^ 2 ~* 1 ~ 1)+ ^ 4 ~* 3 ~ 1 ) + ''' + ^ 2fc ~ i2fc - 1- ^^ 

i 

x • • • (116) 
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Further simplifications are 0:3 = 2 J = 1 and 0:3 — 2 = (—1). As a result the 
expression in (jllfjjl evaluates on the highest weight state to 



2»£(-l 



We also used N(l, 1) = 2 for J = §. The sum over i can be simplified by 
defining variables 

i s = i s — s, s = l,...,2k (117) 
With these variables the expression in ()117f) simplifies to 

2n-2fc is 12 „fc J 

^ A A (n-Jfe)!Jfe! 

*2fc 12=0 ii=0 

In relating the patterns to the symmetrised trace there is a combinatoric 
factor (|3U|). Taking that into account along with the above, we have for spin 
1/2 that the symmetrised power of the quadratic Casimir C n is 1 

C n = ^l±^^ = (2n+1) , f(n) (118) 
( 2n ) ! f^ ((n-A:)!) 2 

Now we can calculate the Lagrangian. Let C = a^ai. Then 



-Sti\ 1+CrH 1-Cs 



-.2 



= _ Str J2 r 4l s 2m (C) l+m (-l) m 

1=0 m=0 ^ ' ^ ' 

= E E r 4l s 2m f(l + m) ( 1/2 ) (^ 2 ) (-l) m 

i=0 m=0 ^ ' ^ ' 

n=0 Z=0 V / V / 

The factor of N(= 2) comes from the trace, and the binomial factors are 
expressed in terms of Gamma functions as 

l/2\ r(3/2) . . 

1 (120) 



/ J r(3/2-/)r(/ + i) 



1 Previous versions of this paper had a missing (2n — 2k)\ on the left resulting in a more 
complicated f(n). The corrected analysis yields simpler energy functions, with the same 
general features. 
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The energy can be computed to be 



E(r,s) = -^f:f:/(/ + m)(2m-l)( 1 / 2 )(^ 2 / ) 



1=0 m=0 

oo n 



= -N^j^f{n){2n-2l-l)( l [ 2 \( 

n=0 1=0 V / V 



1/2^ ( 1/2 

?2 — I 



-l) n ~V 4/ s 2n ~ 2Z 
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If we sum over powers of r first and then over powers of s we get 

Fir «\ - 2 il±^l Ff 1 3 + 4r4 - ( 1 + 2r4 ) ^ f1 oo. 

E(r ' s) " 2 V /[TTH) 2 l( 2' 2T^' (2T^) ' s } (122) 

If we sum over powers of s first and then over powers of r we get 

E ^°^ 2 JT^^ 2 ^T^A- ri) (123) 

These are both useful expressions. 

It is useful to recall the following facts about the hypergeometric function 
2-Fi(z)- This function has branch points at z — 1 and z = oo, and there is a 
branch cut extending from 1 to oo. Formulae are available for the disconti- 
nuity across the branch cut. In ()122j) attempting to continue past s = 1 runs 
into the branch cut. Continuation past s — 1 to s > 1 requires specification 
of the sheet. Physically we are not interested in these super-luminal speeds 
in any case. In ()123|) the hypergeometric function can be continued past 
r = 1 to large positive r even though the the original sum over powers of r 
diverges for r > 1. 

It is instructive to specialize ()122j) or ()123|) to the case of r = or s = 0. 



E(r = 0,s) = 2 (i _ g2)3/2 (124) 

E(r,s = 0) = r-j==4 (125) 
V 1 + r 

From ()125|) we see that the energy grows monotonically with r. This 
means there is no smooth bounce. From ()124|) there is no extremum as a 
function of s at r = 0. Hence there is no exotic bounce at r = 0. Further 
numerical investigation shows that there is no such extremum when r > 0. 
Hence the exotic bounces of the kind studied in section 4 with the first 1/N 
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correction, do not occur for spin 1/2. Whether such extrema occur for other 
finite values of N is an interesting open question. 

The appearance of the limiting velocity s = 1 in (|122|) - (|124j) shows that 
the symmetrised trace prescription correctly captures the expected relativis- 
tic barrier. This does not seem evident from the original form of Lagrangians 
such as dH). 

8 Regimes of validity 

Neglecting higher spatial derivatives implies that the D2-brane picture only 
works for distances larger than the string scale, or r > 1 in dimensionless 
units. We also need to consider higher orders in time derivatives. 

We first discuss issues of regimes which are related to radial velocities 
and their comparison to C, before considering accelerations. For the original 

unperturbed solution, described in Section 2, E = y^jpj and if r is the 

4_ 4 

initial position where r = 0, then f 2 = r °~^ . The subsequent analysis 

of 1/N corrections has shown that higher orders in the 1/N expansion are 
important when 7 2 ~ yC ~ N. This gives a lower bound on r for the large 
iV zero-brane description to remain valid - 

r 4 > | - 1 (126) 

This lower bound can be much less than 1, so we can follow the evolution 
for distances much shorter than the string scale. However it appears that we 
can follow it all the way to zero if r 4 scales like a power of N which is smaller 
than one. 

Indeed, for such choices of r , the paths which solve the equations of 
motion coming from Ci and £2 have a similar behaviour to the leading order 
paths. They describe 2-branes collapsing to zero size with maximal proper 
accelerations which grow as r grows. The maximal speed reached at zero 
radius is very close to the speed of light, with 7 ~ Tq. However when r is 
larger than ~ iV 1//4 , the paths described show qualitative changes. Whether 
we work with C\ or C 2 , we find that there is an upper bound on 7 which is 
much less than r 2 ,. 

8.1 Proper Acceleration and effective mass 

Two useful quantities to study are the proper acceleration and the effective 
squared mass. The proper acceleration a is given by a 2 = Y£^r. Useful 



39 



expressions for a are 



d 2 r 3 ds 3 E r 



E r 



(127) 



This is the relativistically invariant acceleration and should be expected to 
control higher derivative corrections to the Born-Infeld action. For recent 
discussion on the geometry of corrections to brane actions see pJUl 13 1| - An 
effective mass can be defined by m 2 = E 2 —p 2 . For the leading order solution 
one has m = y/1 + r 4 . With this effective mass, the formulae for the energy 
and momentum look like the standard ones of special relativity 
For the leading large N solution 



initial position ro, which is taken to be much larger than 1. It grows to order 
1 at r ~ r 2 / 3 and reaches a maximum of a ~ r 2 , when r = 1. It drops back to 

— 2/3 

order 1 at r ~ r and finally to zero when r = 0. Near the extrema of the 
phase space curves obtained from Ci, where ^ is infinite, the acceleration 
and proper acceleration continuously approach infinity. This suggests that 
the full stringy degrees of freedom are important, or at least we need some 
information about the nature of the terms in the effective action involving 
proper accelerations and higher derivatives. 

For the cases of spin half, and the Lagrangian £ 2 , we can study the 
acceleration numerically. For £ 2) with Uq <~ C 1//4 , the proper acceleration 
behaves similarly to that for the solutions to the zeroth order equations of 
motion coming from £q, reaching large values around r — 1. Somewhat 
surprisingly, for very large Uq, we find that the acceleration remains small 
(less than 1 in magnitude) along the trajectory. This is another indication 
that for such high Uq the 1/N corrected Lagrangians behave very differently 
from the leading order ones. The fact that the proper acceleration remains 
small suggests that there may be interesting time- dependent physics in string 
theory which is non-perturbative in the 1/N expansion but can be reliably 
described by the Born-Infeld actions, nelecting higher derivative corrections 
and massive string modes. 

For the 1/C corrected Lagrangian, Ci, we find that there are regions of 
phase space where the effective mass becomes imaginary. Numerical studies 
show that as we follow the trajectory of a large Z)2-brane down to small 
radius in the space parameterised by (r, s 2 ), we have a minimum radius and 





This shows that the proper acceleration starts at a small value a ~ — at the 
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then the curve continues along increasing s 2 and approaches infinite radius 
as it asymptotes to a finite upper s 2 = s 2 E < 1. The acceleration and proper 
acceleration approach zero as se is approached. The effective mass remains 
real at the extremum but becomes imaginary near se- 
For C\ the effective mass squared is 

m 2 = E 2 -p 2 = U 2 (l + ( 7 2 -C/- 2 ) 2 

-^(7 2 - tT 2 ) 2 (15 7 4 - 16 7 2 + 2 7 2 fT 2 - fT 4 )) (129) 

The region of interest is at large U. Assuming large U and large 7, the 
coefficient of U 2 in the above expression is 1 + ^ — j^yz- This changes sign 
at 7 4 = (7/2. Note that this value of 7 is larger than the large C location 

of the extremum 7 ~ \ consistent with the numerical evidence that the 
extremum is reached before the tachyonic behaviour of the mass squared. 
For fixed initial conditions, there is also a tachyon appearing with £ 2 but 
this also appears at a larger speed than the extremum seen with C±. 

As discussed previously, there are two possible classical evolutions after 
the extremum is reached. One involves bouncing back along the original 
path, with a discontinuity in velocity. The other is to bounce back along the 
branch with increasing s 2 . Here we are finding that the first bounce does 
not encounter the tachyonic region whereas the second does. The immediate 
neighbourhood of the extremum is free from the tachyon but involves infinite 
proper accelerations. The correct string description may require string field 
theory. 

Plots of the effective mass as a function of (r, s) continue to show the radial 
mode becoming tachyonic in certain regions, for £2 as well as for the spin 
half Lagrangian. Following the classical trajectory of the radius as a function 
of 7 for £ 2 , zero radius is reached at finite speed. The effective mass becomes 
imaginary at some point along the classical trajectory if U >~ C 1 / 4 . The 
effective squared mass can be plotted as a function of (r, s) for the spin half 
case. For r = or small, it becomes tachyonic around s = 0.6. 

A summary of the previous comments is that there is a class of pertur- 
bative classical paths which are modified by small 1/(7 corrections. They 
describe the collpase of branes starting at rest at some U which corresponds 
to 7 <~ (7 1 / 4 . In the leading Lagrangian Co the limiting 7 ~ 

C l/4 

are 

reached for Uq ~ C 1 / 4 . The 1/(7 corrections relevant to such paths can be 
computed analytically along the lines of Section (4.3) and (4.4). The exis- 
tence of interesting features such as the extremum we found with C\ requires 
finite iV investigations, which we initiated with the spin half case. If we find 
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such extrema in the finite N case, understanding the complete picture re- 
quires information about higher orders in time derivatives, which should be 
obtained from string theory. This is because the proper acceleration diverges 
at these extrema. Alternatively we may go beyond the framework of effective 
actions and study the problem in string field theory. It may also be inter- 
esting to explore what happens to such extrema when one considers simple 
higher derivative Lagrangians which are designed to put an upper limit on 
the proper acceleration [32~] . 

8.2 The physics of the tachyon 

A heuristic argument can be used to predict the existence of the tachyonic 
effective mass which we found. Imagine an observer sitting at the north 
pole of the collapsing spherical D2-brane, at $ 3 = R. In the rest frame 
of this observer, the north-south axis is Lorentz-contracted, so that the 2- 
brane looks more like the surface of a pancake. The spherical brane has 
zero net D2-charge, and can be thought as having positive charge at the 
north pole and negative charge at the south pole. This follows since the 
D2-charge density is proportional to [$i,$2] ~ ^3, which is positive at 
the north pole and negative at the south pole. Locally then, our observer 
at the north pole sees themselves at rest on a D2-brane with an anti-D2 
brane approaching at high speed. We know that the D2-anti-D2 system has 
a tachyon when the separation is close to the string scale. The observer 
also sees a density of magnetic flux which is a distribution of DO-charge, but 
D0-D2 systems also have a tachyon. It would be interesting to make a more 
quantitative connection between this picture and the tachyonic effective mass 
obtained from C\, C2 and the spin half case, in order to understand better 
this "pancake tachyon" . The insights from E3J EHj which discuss brane- 
anti-brane systems in the presence of motion and flux may be useful. Another 
situation where tachyonic behaviour of a radial brane variable has been found 
recently is described in 

It might be argued that the tachyon is an artefact of the symmetrised 
trace prescription, which is not the correct supersymmetric non-abelian DBI 
theory. It has been shown that the symmetrised trace prescription must 
be corrected, based on BPS energy formulas |371 EB1 EH| • However, we find 
it likely that these tachyonic features will survive with these corrections. 
Indeed it can be argued that L\ is not modified. The form Str(aa) n = 
N(C n +a 1 (n-l)C n - 1 + n(n-l)(n-2)(an+b)C n - 2 + . ..) follows from general 
arguments which will hold true even when the symmetrised trace prescrip- 
tion is replaced by something more complicated such as weighting different 
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symmetry patterns with different coefficients. For example, we know that 
the first correction of the form C*( n_1 ) must vanish for n = and 1 because 
for these values of n we can simply evaluate the traces and check that the 
leading term C n is accurate. For n — 1, we have tr^aia/) = NC . Assum- 
ing that the n-dependent coefficients are sufficiently nice functions of n, it 
would follow that the first correction is unchanged by these corrections. The 
failure of the symmetrised trace prescription starts at tr(F 6 ) which trans- 
lates here into tr(ctiai) 3 . Once the structure of the tr(F 8 ) terms are known 
the numbers a, b appearing in the coeffient of C n ~ 2 can be determined, and 
will likely be different from the values 7,-1 of (J23*j) . We have checked that 
the effective mass to order 1/C 2 can become imaginary at large 7 for any 
choice of a, b. It is still possible to argue that the tachyon seen in the 1/N 
expansion is really an artefact of the failure of the 1/N expansion itself, and 
that the correct supersymmetric non-abelian Born-Infeld at finite N will not 
allow such tachyonic behaviour. An important future direction is to develop 
a proof (which does not look at all obvious) along these lines which disposes 
of the tachyon, or to develop a more concrete quantitative formulation of the 
physical nature of the tachyon as outlined above. The expectation that the 
1/N expansion captures some qualitative features of the finite N physics sug- 
gests that the latter avenue will be more fruitful. We have also seen that the 
symmetrised trace prescription does lead to formulae for the energy which 
correctly capture the expected relativistic limit s — 1 in the spin half case. 



9 Summary and outlook 

A large iV approximation to the non-abelian zero-brane action of Myers for 
time-dependent fuzzy sphere configurations gives equations for the radius 
which agree with the expected dual picture in terms of a D2-brane with a 
magnetic flux. These equations have the same form as those for a relativistic 
particle with a position dependent mass. The 1/N corrections to the zero- 
brane action give rise to modifications of these equations. We studied the 
energy function for these modified equations in detail. 

The simplest solution of interest at large N, with Lagrangian Co, describes 
a D2-brane, with initial radius R ^> l s collapsing to zero radius. Classically 
this solution can be patched with an expanding solution at zero radius, but 
the velocity is discontinuous at that point. The zero brane description allows 
us to deduce that the collapsing radius can be trusted to distances less than 
the string length. However the speed at very small radius for such a large 
initial membrane is close to the speed of light. At these speeds, higher orders 
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in the 1/N expansion of the zero brane action become important. Another 
feature to bear in mind is that the proper acceleration becomes large along the 
classical trajectory. This raises an interesting question of whether the form 
of the time-dependent solution can be argued to be unchanged throughout 
the trajectory from large R down to zero, perhaps using arguments along the 
lines of 0Dj. 

Using the first order in 1/N corrected Lagrangian L\ , we showed that for 
a large range of initial conditions R > l s (r > 1 in dimensionless variables), 
the classical path obtained from Ci is qualitatively similar to that obtained 
from C . However for initial conditions which allow 7 to reach near 

£1/4 

along the path, the path described in phase space encounters a minimum 
at non-zero speed. This conclusion can be reached from the analysis of 
the energy as a function of r and s, where r is the dimensionless radius 
and s the dimensionless velocity. The contours of constant energy starting 
from zero velocity at some large radius cannot be continued to zero radius, 
but rather have a minimum. At the minimum, = 0, while s m i n 7^ 0. 
Such minima do not occur in simple mechanical systems where E = + 
V(x). This occurs in our equations because of the peculiar mixing of the 
coordinate and velocity in the energy function. At the extremum, the velocity 
is close to the speed of light. As N — > 00 this extremum recedes to the 
relativistic barrier s = 1, which is why we do not see it in the leading large 
iV approximation. The bounce at the extremum is somewhat peculiar, it 
involves a discontinuity in the velocity and also an ambiguity in the choice 
of trajectory after the bounce. These features are made clear by a local 
analysis of the energy function in the neighbourhood of the singularity. We 
discussed the quantum mechanical behaviour of the system near this local 
singularity, and argued that quantum mechanics would provide probabilities 
for the two trajectories. The discussion is not complete since the equations 
we are trying to quantise have a non-linear dissipative nature, hinting that 
the correct quantum treatment requires extra degrees of freedom, such as 
those that naturally occur in the string set-up when we go beyond the low 
energy effective action approach and include higher string excitations. 

Another consequence of the fact that the velocity at the extremum is 
close to the speed of light is that the 1/N correction term to the Lagrangian 
is comparable to the leading term for such velocities. This means that the 
1/N approximation is breaking down in the neighborhood of the extremum. 
This motivated an analysis of higher orders in the 1/N expansion. At the 
second order in the 1/N expansion we found that the extremum disappears, 
while it reappears again at the third order. 

This led us to study the exact Lagrangian for finite N. In the case N = 2, 
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we found that the extremum does not exist, but for general N the existence 
or otherwise of such extrema is a very interesting open question, which could 
be addressed using the techniques developed in this paper. If such extrema 
exist then we will have a new mechanism in string theory for a brane bounce 
which is intimately related to the non-Abelian structure of D-brane actions. 

The technical computations of the 1/N corrections require care in dealing 
with the symmetrised trace. This symmetrised trace leads to the evaluation 
of certain £0(3) invariants. We presented two approaches to the calculation. 
One proceeded by evaluating on highest weights. The other used a diagram- 
matic approach for re-ordering the various terms in the symmetrised product. 
This latter approach makes interesting connections with knot theory, which 
may provide a source of new techniques and ideas in this context. 

Further computations with C\ included a perturbative computation of 
the time taken to collapse from an initial radius to some final radius. This 
perturbative approach is best suited to initial and final radii which are not 
too far apart. We also computed the energy as a function of momentum 
giving an interesting deformation of the dispersion relation for a massive rel- 
ativistic particle, derived directly from the non-Abelian structure of the zero 
brane action in string theory. Whereas many string-inspired deformations 
of relativistic dispersion relations have been discussed, this deformation is 
the first that follows directly from the non-abelian nature of stringy D-brane 
actions. 

In Section 7, we gave a detailed discussion of regimes of validity for the 
different approximations, taking care to use the proper acceleration rather 
than the ordinary acceleration as a measure of whether higher derivative 
corrections are important. Quite generically, we found that the generalised 
effective mass defined by a/ E 2 — p 2 became tachyonic in certain regions of 
phase space close to the speed of light. We discussed the physics of this 
"pancake tachyon" as associated with a geometry which looks very similar 
to system of brane and anti-brane. We gave arguments against attributing 
the tachyon to the inadequacy of the symmetrised trace prescription used in 
the non-abelian Born-Infeld. 

Here we have considered the simple situation of a spherical D2-brane 
collapsing from a dual DO-point of view using the fuzzy 2-sphere construction. 
It will be interesting to explore similar brane collapse phenomena for higher 
branes using higher dimensional fuzzy spheres [UJ H21 HH 113 ESI HZj- 
This may be a new way to explore cosmological bounces in a braneworld 
context developing works such as [48J, with the additional ingredient of the 
non-abelian Born-Infeld action. Position dependent effective masses have 
been considered in |49| 1571] and earlier references therein. The leading order 
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Lagrangian Co is a simple example of a relativistic system with position 
dependent mass. An effective mass can be defined using m 2 = E 2 —p 2 for the 
higher order Lagrangians. We found tachyonic behaviour in certain regions 
of phase space close to the speed of light. This problem of the collapsing 
brane may also have similarities with gravitational collapse of thin shells jSI] . 
Another direction is to look for a gravitational background which corresponds 
to the time-dependent system of large N zero branes, and find a spacetime 
interpretation for some of the features we have found such as the bounces 
and the tachyonic effective masses. 
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